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We study single-electron transport in a magnetic molecule by using a stationary rate equation 
approach. In the molecule, two magnetic ions and are connected to each other through a third, 
nonmagnetic ion. The magnetic ions are coupled to ideal metallic leads and a back gate voltage is 
applied to the molecule, forming a field-effect transistor. A Hamiltonian which includes inter-ion 
hopping, on-site repulsion, and magnetic anisotropics is employed to describe this molecule, resulting 
in an energy spectrum similar to that of single molecule magnets in the giant-spin approximation. 
An external in-plane magnetic field is then used to drive the molecule to a diabolical point, where 
states with maximum total spin with opposite directions are degenerated. Both linear and nonlinear 
transport are evaluated near the diabolical point, showing features that can be attributed to Berry- 
phase interference of spin tunneling paths. 



I. INTRODUCTION 

From the various types of magnetic states in matter, 
ferromagnetism and its microscopic causes is one of the 
most intriguing topics. A very singular class of mag- 
netic systems are single-molecule magnets (SMMs)p3ES 
Molecules in this class typically have a large net spin 
ground state and exhibit unusual attributes su ch a s quan- 
tum tunneling of the mag netization (QTM)pEl a rela- 
tively large decoherence timeP and Berry-phas e in terfer- 
ence effects in the presence of magnetic fields!^ Com- 
monly, SMMs are composed of transition metal ions 
(open 3d or 4/ shells) bridged by ligand atoms and 
molecules. The Hamiltonian that describes well the mag- 
netic properties of SMMs is the so-called giant spin ap- 
proximation (GSA) model. For instance, in its simplest 
form, the GSA Hamiltonian can be written as 

H GSA = -DS 2 Z +E (S 2 x -S 2 y ), (1) 

The net magnetization has a preferential direction (an 
easy axis z in this example) as a consequence of the 
dominant uniaxial anisotropy (D <C E). The in-plane 
transverse anisotropy allows the total spin to transit be- 
tween different values of S z . In this particular example, 
the molecule has a predominant second-order anisotropy 
signaling a rhombic symmetry. 

A more microscopic description of SMMs comes from 
considering interactions at the ion level. In this case, one 
uses instead a multi-spin Hamiltonian of the type 

H multi-spin — ^ [ d{S z + &i [S x ^y^)1 ^ JijSi • Sj, 

i ijtj 

(2) 

Here, each magnetic site of the molecule has itself local 
uniaxial and transverse anisotropies di and e,, respec- 
tively. In addition, there is an effective ferromagnetic 
interaction between pairs of sites parameterized by Jij. 
Both models have been compared in Ref. |9] For zero 
transverse anisotropies (ej = E = 0), the diagonaliza- 
tion of the Hamiltonians in Eqs. ([!]) and ^ yield a 



two-fold degenerate ground state involving spin states 
parallel (S z — S) and anti-parallel (S z — —S) to the uni- 
axial anisotropy axis. If the transverse anisotropy terms 
(ej and E) are nonzero, then rotational S z symmetry is 
broken and the two-fold degeneracy is lifted, with the 
ground and first excited states being antisymmetric and 
symmetric combinations of the S z = ±S states. In the- 
oretical studies, due to the large Hilbert space involved 
and the strong interplay between the many microscopic 
parameters, it is very challenging to take into account 
explicitly the overlap between orbitals and their sym- 
metries in magnetic molecules. One often employs the 
empirical Goodenough-Kanamori rule d that dictate 
the nature of interaction between magnetic ions. In addi- 
tion, one also includes an effective interaction intermedi- 
ated by diamagnetic atoms (the Anderson superexchange 
interaction) P21 

In recent years, interest emerged in exploring how the 
magnetic properties of SMMs may affect the molecule's 
electronic transport properties. Effec ts due to the 
QTM^and Berry-phase interf erenc e ! 14 * 15 ^ and their inter- 
play with the Kondo resonanc e! 16 * 17 ! have been proposed. 
Several authors have also used first-principles calcula- 
tions to investigate the electronic configuration, the mag- 
netic properties, and coherent transport in SMMsP^HZU 

In this paper, we investigate the dependence of sequen- 
tial tunneling transport on an applied transverse mag- 
netic field in a three-ion model that contains the essential 
microscopic details necessary for reproducing a SMM be- 
havior. The microscopic model comprises two magnetic 
ions bridged by a third diamagnetic ion. The microscopic 
model takes into account the valence, ligand fields, and 
orbital energies of the ions, as well as direct and exchange 
Coulomb interactions present in the molecule. Our goal 
is to develop a simple yet phcnomcnologically complete 
microscopic model of a SMM that goes beyond the gi- 
ant spin Hamiltonian model and is amenable to realistic 
electronic transport studies. 

We explore the energy spectrum of this system for sev- 
eral points in the parameter space of the Hamiltonian. 
We find that for a certain parameter range a magnetic 



2 



bistability, similar to that of SMMs, develops. We then 
focus our attention on Berry-phase interference effects 
on the single-electron transport through the molecule 
since this is one of the most unusual features expected 
of a SMM. In our calculations, transport is evaluated 
under an applied transverse magnetic field whose magni- 
tude and direction are tuned with respect to a diabolical 
point (see Ref. 22). At the diabolical point, Berry-phase 
destructive interference of spin paths of the molecule is 
maximal, completely compensating the mixing of oppo- 
site spin orientations brought by the in-plane anisotropy. 
As a result, the ground-state degeneracy of spin states is 
restored. . 

The paper is organized as follows: In Sec. [TTJ we de- 
scribe the model Hamiltonian for the three-ion magnetic 
molecule. In Sec. |III| presents our choice of model pa- 
rameters and the energy spectrum and symmetry of the 
states of the molecule in the absence and in the presence 
of transverse anisotropy and magnetic field terms. In Sec. 
|IV| we evaluate the incoherent sequential current and dif- 
ferential conductance of the molecule for zero field and 
for a nonzero field tuned near the diabolical point. In 
Sec. |VJ we discuss the effects of Berry phase interference 
on transport, which are most visible in the molecule's 
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differential conductance. Finally, in Sec. VI we analyze 
our results and draw some conclusions. 



II. THE THREE-ION MODEL OF A SMM 

A minimal realistic molecular core capable of repro- 
ducing the main features of a SMM such as large to- 
tal spin (S > 1/2), uniaxial and in-plane anisotropies, 
and Berry-phase interference consists of two transition 
metal ions bridged by a diamagnetic ligand ion. Con- 
sider the system shown schematically in Fig. [I] The 
two transition metal ions a and c have a 3d 8 electronic 
configuration and a total spin S = 1 each. They in- 
teract magnetically through a superexchange interaction 
intermediated by a diamagnetic O 2- ion (electronic con- 
figuration 2p 6 ), represented by b in Fig. 1 The five-fold 
degeneracy of the 3d orbitals in the magnetic ions is bro- 
ken due to the bonding to ligands. To simulate such an 
effect, we assume that a weak orthorhombic field acts 
on a and c, inducing identical local uniaxial and trans- 
verse anisotropies on each site. Thus, for an orthorhom- 
bic symmetry (I)ih point group), the ground states of 
ions a and c have two occupied unpaired single-particle 
orbitals (eii % and c± 2) respectively) which are symmetric 
and antisymmetric combinations ofd x 2_ y 2 and d 3z 2_ r 2 d 
states. We also consider a 90° angle between their bonds. 
As a result, the d x 2_ y 2 components of the single-particle 
orbitals in both ions a and c overlap separately with the 
p x .y orbitals &1.2 of ion b. Under this assumption, we use 
the appropriate Slater-Koster two-center integrals, which 
then yields E Pxidx2 _ y2 = E Py ^ 2y2 = (pdcr)v / 3/2. This 
configuration favors a ferromagnetic superexchange in- 
teraction. In our case the a two-center integrals are zero. 




Figure 1: (Color online) Scheme of a simple three-ion SMM. 
The magnetic ions labeled by a and c interact through the 
diamagnetic ligand ion b with a bond angle of 9 — 7r/2, so that 
the d x 2_ y 2 orbitals (green/solid lobes) overlap with thep x and 
p y orbitals (purple/dasher lobes) of the ligand separately. 




Figure 2: Scheme of the electron hopping (double-arrow 
lines) between magnetic ions and the ligand's orbitals. For 
a 3d 8 ion in a weak orthorhombic ligand field, the unpaired 
single-particle orbitals (say for ion a) are ip a2 = Kid x 2_ y 2 + 
K2rf 3z 2_ r 2 and tpa-i = ki d 3z 2_ r 2 — k,2 d x 2_ y 2 , where Ki + K2 = 
1. We regard K2 as a small parameter: \k2\ *C 



The scheme of these selective orbital overlaps is shown 
in Fig. [2j The features of our model are captured by the 
effective Hamiltonian 

-^molecule = H a + Hb + H c + H a b + Hb c (3) 

where 

i=l,2 i=l,2 

F U a ^ ^ai ,(7 ^ft2 d a S ai • S a2 

(7,(T'=t,4. 

- d a S z>a + e a (Sl a — Sy a ) , (4) 
with a = a,b,c and (4 = = 0, 

Hab = t ^2 (4 2 ,a c 6i,^ + H.c.) 
CT =t4 

+ *' E (ci^Q^+H.c), (5) 



and 

H bc = t (cL^.+H.c.) 

<7=t4 

+ < ' E ( c L, ff c Cl , ff '+H.c). (6) 

In Eqs. (j3]), ([HI, and ([H]), n ai)<7 = c^.^Cq,^, where cJ, i>(T 
(cctjo-) creates (annihilates) an electron with 2 spin pro- 
jection a in the orbital i of ion a and follow the standard 
fermionic anticommutation relations. The total spin op- 
erator associated to the orbital i in ion a is S ai while 

In Eqs. ([5j| and (g, we include both spin preserving 
and spin flipping hopping terms. In Eq. Gj), the first 
term accounts for the filling of the single-particle orbitals, 
while the second and third terms represent the on-site 
Coulomb repulsion. The fourth term enforces Hund's 
first rule, maximizing the total spin on each ion (J a > 0), 
and the last two terms are the uniaxial and transverse 
on-site anisotropies produced by the rhombic ligand field 
environment and by the spin-orbit coupling within each 
magnetic ion. Typically, \d a \ > \e a \. 

III. ADJUSTING THE MODEL PARAMETERS 

The model Hamiltonian in Eq. ^ has a large number 
of parameters that have to be properly adjusted in or- 
der to produce the phenomenology expected of a SMM. 
Below, we discuss our choices of parameter values. We 
first consider the Hamiltonian in the absence of trans- 
verse anisotropy (e a = and t' = 0). Then we include 
transverse anisotropy and the coupling to an in-plane ex- 
ternal magnetic field in order to produce diabolical points 
in the spectrum of the molecule. 

A. Total S z spin in the absence of transverse 
anisotropy and of spin-flip hopping 

The eigenvalues and eigenstates of the Hamiltonian in 
Eq. ([3]) are very sensitive to the values of its parame- 
ters. A convenient way to produce a spectrum similar to 
those of SMMs is to start with S z as a good quantum 
number by setting e a = t' = 0. Then, one considers how 
the eigenstate amplitudes resemble the Clebsch-Gordan 
amplitudes obtained for two 5 = 1 particles when dif- 
ferent choices of parameters t, s ai , U a , U' a , J a , and d a 
are used. Due to the molecule's symmetry, we can as- 
sume e ai = e Ci = e i} U a = U c = U, U' a = U' c = V, 
Ja — Jc = J, arL d d a = d c = d. For simplicity we have 
set Ub = U' h = 0. The single-particle orbital energies of 
the bridging ion are set as e&j = Eb 2 = £&. 

Realistic values for some parameters can be obtained 
from the review by Imada, Fujimori, and TokurgP^l a nd 
the recent work by Kim and MinPS We set the on-site 
and intraorbital Coulomb interactions to U — 7 eV and 
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U' = 4.5 eV, respectively. The Hund's rule parameters 
are J = 0.95 eV and Jb = 1.3 eV. The orbital energies 
have the values e\ = 0.71 eV, e 2 = 0.72 eV, and = 0.05 
eV. 

The uniaxial anisotropy parameter is set to d = 58 fieV 
(0.67 K). For this particular choice of parameter values, 
we calculate the total S z spin of the molecule. Figure 
[3] shows \S Z \ for the ground state (two-fold degenerate), 
first excited state (two- fold degenerate), and second ex- 
cited state (nondegenerate) versus the hopping parame- 
ter t. It is clear that for the range of values probed, the 
ground state has the highest possible total spin, 5 = 2. 
For t > 1.2 eV, the states with total spin 5=1 be- 
come the first excited states, while the state with total 
spin zero moves up to the second excited state position. 
This is the typical case for a SMM within the description 
provided by the giant spin Hamiltonian of Eq. ([!]) . 
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Figure 3: (Color online) Total S z spin of the molecule. For 
< t < 1.19 eV, the lowest energy level is a ± 2} state, 
the first excited state has S z = 0, while the next (degener- 
ate) excited states have S z = ±1. We find the typical SMM 
behavior in the 1.2 eV < t < 1.5 eV range, where the first 
excited states have S z = ±1 and the highest energy level has 
S z — 0. The crossing in the graph denotes the point where 
total S z spin changes for the first and second excited eigen- 
states. For t = 1.5 eV, the energy splitting between the two 
lowest eigenstates (S z = ±2 and S z = ±1 ) is of the order 
of the uniaxial anisotropy parameters d a , with the value of 
A = 57/ieV, while the splitting between S z = ±1 and S z = 
is A = 8.43 ^eV. 



B. Anisotropy and diabolical points 

We now consider the inclusion of anisotropy terms, 
which make the eigenstates to be symmetric and anti- 
symmetric combinations of \S Z ) states. In addition, wc 
consider the effect of an applied transverse (in-planc) 
magnetic field on the molecule's energy spectrum. The 
transverse field Hamiltonian is given by 

-ff field = b I 

COS (p + 9y,a Sy^ a Sill (j) \ , (7) 
\a—a,c / 
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where b = hb\B\ and B = (B cos <j>, B sin (j>, 0) is the 
external magnetic field. Following Refs. [25] and [26j a 
reasonable estimate of the anisotropic ^-factors would be 
9x,a = 2.270 and g y a — 2.269. We set the spin con- 
serving hopping parameter to t = 1.5 cV, while the spin- 
flipping and transverse anisotropy parameters are set to 
t' = 9/ieV and e a = e c = 0.17d, respectively. We do 
not consider coupling of the magnetic field to the middle 
ion since it would change the energy levels of the high 
energy sector by a very small amount, given that b <C U. 
In Fig. [4j we show the resulting splitting of the two low- 
est energy levels versus the angle 4>, while in Fig. [5] the 
splitting is shown as a function of the applied field magni- 
tude b. The angle tj> = tt creates symmetric spin paths for 
the tunneling of the molecule's magnetization, leading to 
maximum destructive interference. It is at this particular 
direction that we plot the energy splitting in Fig. [5] The 
periodic modulation of the splitting with b is a character- 
istic feature of SMMs (see Refs. 14171270 . The separation 
between successive zeros of the splitting for a system de- 
scribed by the GSA [Eq. Q] under a field /ifi e id = —Qx^ 
(with b = hbB) is given by A& gsa = {2/g x )^2E{E + D). 
Note that using this expression and our choice of values 
for the anisotropy parameters e and d, the predicted pe- 
riod of the oscillation would be A6 gsa = 32.2 fieV, which 
is not in agreement with the actual value of Ab ps 8.5 /LteV 
obtained from Fig. [5j The derivation of an analytical 
expression for the period of the splitting may not be pos- 
sible for our Hamiltonian and we do not attempt it here. 
In the following, we will the fix the direction of the field 
to <j) — 7r and only consider the dependence with respect 
to the field magnitude b. 

We now consider how the symmetry and energy of the 
eigenstates are affected by changes in the magnitude of 
the magnetic field. For this purpose, we define the sym- 
metry coefficient 
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C 2 



(8) 



where are the amplitudes of the two degenerate 

states along the S z = ±2 basis states. If Cs = the 
eigenstate is antisymmetric, whereas if Cs = C\ 2 the 
eigenstate is symmetric. The results are presented in 
Fig. [6] We observe that the ground state changes from 
symmetric to antisymmetric at the first diabolical point 
(left dashed line), while the first excited state does the 
reverse. The energy splitting between the linear combi- 
nation of S z states (symmetric and antisymmetric) van- 
ishes at the diabolical points, since the two lowest energy 
states become degenerate at those points. 



IV. SEQUENTIAL TRANSPORT THROUGH 
THE MOLECULE 

We study electronic transport in our model system by 
connecting the molecule to two reservoirs of noninteract- 
ing electrons. The current through the molecule can be 
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Figure 4: (Color online) Splitting of the two lowest energy 
levels versus the transverse field angle <j>. The angle is mea- 
sured with respect to the positive x axis. We see that for 
various values of the magnetic field, the lowest splitting oc- 
curs at <A = 7r. 
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Figure 5: (Color online) Splitting of the two lowest energy 
levels versus the magnitude of the magnetic field b at an angle 
of 4> — 7T. We observe two diabolical points where the splitting 
between the two lowest energy levels goes to zero. The first 
diabolical point occurs at b ~ 8.5 fieV (B w 0.14 tesla), while 
the second one occurs at b « 17 /xeV (B w 0.28 tesla). 



controlled by applying a voltage difference between the 
reservoirs, as well as by changing the total charge of the 
molecule through an applied backgate voltage. In order 
to evaluate single-electron transport, it is necessary to 
diagonalize the Hamiltonian for configurations with dif- 
ferent total charge (i.e., different number of valence elec- 
trons in the ions). For simplicity, we only consider two 
charge states and we define (N = 0, S = 2) as our initial 
state. In this state, the molecule has a 3d 8 — 2p 6 — 3d 8 
electronic configuration. Here, S denotes the total spin 
of the molecule. We then allow one electron to hop from 
the reservoirs into the molecule and restrict it to be local- 
ized either on the a\ or c\ orbitals, bringing the molecule 
to the (N — 1,S — 3/2) state, which comprises both 
3d 9 — 2p 6 — 3d 8 and 3d 8 — 2p 6 — 3d 9 electronic configu- 
rations. We assume that the coupling to the reservoirs 
is such that there is equal probability for an electron to 
land in either one of the two magnetic ions. Since one of 
the ions changes its oxidation state when an electron is 
added, its uniaxial anisotropy terms and the g-factors 
change as well. We take into account some of these 
changes in the Hamiltonian for the S = 3/2 configura- 
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Figure 6: (Color online) Symmetry of the ground state (filled 
squares) and first excited state (empty squares) versus the 
magnetic field magnitude b. At b — the ground and first ex- 
cited states are jgs) = 0.68(|+2) + |-2))+0.27x (contam.) and 
|ei> = -0.69(| + 2) - | - 2>) + 0.22 x (contam.), respectively, 
where (contam.) represents a contribution from spin states 
other than | ±2). Notice that as the transverse magnetic field 
grows, the amplitude Ce tends to decrease, allowing for an 
increasing admixture with S z = states. At the second dia- 
bolical point (right dashed line), the symmetry of the states 
changes again. 



tion. Following Rcfs. [T2]and[2Hl we assume a reduction in 
the uniaxial anisotropy and choose d N=1 = d/2. Notice 
that although the parameter e N=1 is nonzero, it is not 
relevant to the ion receiving one electron since its total 
spin becomes S a = 1/2. For S a = 1/2, the local ground 
state of the ion is a Kramers' doublet, thus e N=1 does not 
take part in the Hamiltonian for the (N = 1, S = 3/2) 
configuration. In the literature,^ there are indications 
that the anisotropic g factors may take a different value 
for the (N =1,5 = 3/2) configuration. However, we do 
not expect the change to be significant and therefore, for 
the sake of simplicity, we assume the same g factors used 
for the (N = 0, S = 2) configuration. 

The general recipe for evaluating sequential, incoherent 
transport is to solve a set of coupled differential equations 
for the time evolution of the quantum state probabilities 
of the molecule. The standard rate equation for time 
evolution of the probability of an arbitrary state is given 



dp a 

dt 



-p, 



7 



7— >a 1 



(9) 



where r Q ,_ >7 is the transition rate between different 
charge states a and 7 and vice versa. The meaning of the 
right-hand side of the equation is the following. The first 
term comprises outgoing terms describing the transition 
from the initial state a to a state 7 by means of taking 
away or adding an extra electron. The second term on 
the right-hand side of Eq. Q represents processes where 
an initial state 7 transitions to the final state a (again, 
by means of taking away or adding an extra electron). 
Equation ^ can be written in a matrix form 

dP 

dt 



TP, 



(10) 



where P is a vector containing the probabilities {p a } and 
the matrix T contains elements {r a7 } such that 



r, 



r^—^a ff 
— J]] r a _> 7 . 

7^a 



a + 7 



(11) 



In the stationary regime, we look for the steady-state 
probabilities p a , this is, when dp a /dt = 0. In matrix 
form, this is achieved by determining the eigenvector P of 
the matrix T with a zero eigenvalue. The transition rates 
from a state a to a state 7 when adding (or subtracting) 
electrons (e~) to (from) the molecule are given by 



ri 



lL n fr(^ ai - eV T ), e in, 

7<L 7 [1 - Ufacr, - eV r)} > e ~ ° ut > 



(12) 



where the index r = L,R refers to the left and right 



reservoirs, respectively. Here, f T {x) = 1/(1 



) is the 



Fermi distribution of the left (L) or right (R) reservoirs. 
The electrochemical potential of the molecule, i.e., the 
energy required to go from a state a to state a 7, is 
defined as 



(13) 



where e a and e 7 are the energy eigenvalues corresponding 
to the states a and 7, respectively, and V g is the backgate 
voltage (we set the lever arm coefficient rj = 1). We 
denote the splitting between eigenstates as A Q7 = e a — 
e 7 . The coefficients 7^ 7 are the tunneling rates between 
eigenstates of the molecule and are given in the Golden 
rule approximation by 



Tq,7 



m, 
out, 



where p is the reservoir's density of states and 



t: 



jo- 



l7> 



(14) 



(15) 



denote the tunneling matrix elements. The operator cj CT 
creates an electron with spin a on the single-particle or- 
bital j. To simplify the calculations we assume that tj jT is 
independent of the orbital, so we drop the j index. Since 
different charge configurations of the molecule involve 
states with different spin quantum numbers, selection 
rules for spin transitions become important in the trans- 
port calculations. The tunneling matrix elements satisfy 
the selection rules \S a — S^\ = 1/2, \S za — S Z1 \ = 1/2 
and \N a — JV~| = 1, for any two eigenstates of the Hamil- 
tonian. The current through the left lead is given by 



Ir 



(16) 



where a, 7 are indexes for different charge states. The ex- 
pression for the current is multiplied by 1 if the molecule 
is receiving an electron and by — 1 if an electron is going 
out to the left reservoir. 



6 



V. BERRY-PHASE INTERFERENCE EFFECTS 
ON THE TRANSPORT 

In order to see signatures of Berry-phase interference 
in the molecule's electronic transport, we consider only 
transitions between the two lowest energy levels of each 
spin-charge sector, as shown in Fig. [7J Since the energy 
levels for the (N = 0,5 = 1) and (N = 1,S = 3/2) 
configurations are linear combinations of \S Z = ±2) and 
\S Z = ±3/2) states, respectively, the corresponding selec- 
tion rules are trivially fulfilled and transitions between all 
states are possible at 6 = 0. 
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Figure 7: (Color online) The two lowest energy levels for the 
N — and N = 1 charge sectors, with total spins S = 2 and 
S — 3/2, respectively. The arrows denote allowed transitions 
for 6 = 0. 

It is worth noticing that although we choose only the 
two lowest energy levels, interesting physics beyond the 
reach of the GSA can be found if all combinations of 
S z states corresponding to a particular total spin S of 
a molecule are considered. This is because anisotropy 
contributions tend to contaminate the lowest energy lev- 
els with a small admixture of S z — states, opening 
transitions in the transport that would otherwise have 
been forbidden by selection rules. An interesting study of 
these contributions to transport can be found in Ref. 1161 
Since we are primarily interested in how the modulation 
of the transverse anisotropy splitting affects transport, 
we only consider the two lowest energy S z states. For 
other states, the energy cost to access different transi- 
tions would be too large compared to the splittings gen- 
erated by the transverse anisotropy. 

Because the splitting of the two lowest states due to 
transverse anisotropy for each charge sector is very small, 
of the order of /xeV (see Fig. [8]) , a very small tempera- 
ture is required to resolve features in the electronic trans- 
port that can be associated to Berry-phase interference. 
Thus, we set the temperature in the reservoirs to T = 0.5 
mK. In order to bring the Berry-phase interference ef- 
fects to higher temperature ranges, we attempted to in- 
crease the tunneling splitting rates by moving slightly 
away from the critical angle where the diabolical point 
develops. Unfortunately, we found that the splittings for 
the two spin-charge sectors grow at the same rate with 
the applied magnetic field, and their difference (which is 
the primary scale controlling the Berry phase interfer- 
ence effects), did not increase significantly. Thus, from 
an experimental perspective, it may be rather difficult to 
resolve all transport resonances in molecules with transi- 
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Figure 8: Splittings between the two lowest energy levels of 
the S = 3/2 (solid line) and S — 2 (dashed line) total spin 
sectors. We observe that As =3 / 2 < As = 2. The first diabolical 
point for S = 3/2 is at 63/2 ~ 1.21 peV (equivalent to 0.02 
T). 



tion metal ion cores. 

For the transport calculation, we tune the gate voltage 



to V g (b) 



M=l 



(6) so that the ground states 



energies for both spin-charge sectors are aligned. We set 
the density of states of the reservoirs and the tunneling 
coupling to arbitrary constant values. Transport is eval- 
uated at b — and b = 6 3 / 2 , the latter denoting the 
magnitude b$ of the field for which a diabolical point oc- 
curs for the spin-charge sector with total spin S = 3/2. 
The left-lead current 1^ and the corresponding differen- 
tial conductance (dlL/dV) are shown in Figs. [9] and 10 
respectively (the subscript L is omitted in the plots). 

When the magnetic field is zero (b = 0), transi- 
tions between excited states and between ground to ex- 
cited states in different sectors can be appreciable in 
the last three steps of the current in Fig(9j which oc- 
cur at V b w 33.5,45.8,79.9 (in units of k B T/e). The 
steps in Fig. |9| correspond to the linear conductance 
peaks in Fig 
S 
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As the magnetic field is tuned to the 
3/2 diabolical point, the degeneracy of the S z = 3/2 



states is restored, suppressing the | e^ v= ) O | = ) and 
|gs 7V=0 ) o |e^ =1 ) transitions. This suppression is clearly 
visible in gray plots of dl/dV versus bias and gate volt- 
ages (Figs. [Tlj and [l~2"| . In Fig 11 the positive and neg- 
ative slope lines corresponding to the |ej V=0 ) <-> 



and |gs ) -H- |ej ) transitions, respectively, are elim- 
inated, as shown in Fig. [12] The effect of the field there- 
fore, is to allow current to go in or out of the molecule 
only if it is between ground states of the molecule or if 
I gs 7V=1 ) O I e^ v=0 ) takes place. The latter corresponds to 
the surviving positive slope line in Fig. |12| 

In Fig. [TUJ we also notice that a shift in the con- 
ductance peaks occurs for the transition between ground 
and first excited states of the N = 1 and N = charge 
states. This is because the splitting of the N = charge 
sector also varies with the application of the field (see 
Fig. [8]) . The field also causes a shift of the charge degen- 
eracy point in Fig. |12| Since our choice of temperature 
is rather small, the increments of the gate voltage in the 
gray plots are very small. Thus, to improve visibility of 
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Figure 9: Current through left lead versus bias voltage at 
6 = (solid line) and 6 = 6 3 / 2 (dashed line) for T = 0.5 
mK (current shown in arbitrary units). The bias voltage is 
plotted in units of ksT/\e\. We observe that the resonances 
at Vt = ±33.5, ±45.8 (in units of fc_gT/|e|), which are present 
for 6 = 0, are eliminated due to the Berry-phase interference 
when b = bg/2 



the data, we shifted the gate voltage. The numerical val- 
ues for the gate voltage axis correspond to the rescaling 
e(V s - V s )/k B T where eV s = 13.175 eV. 



VI. CONCLUSION 

We have studied the incoherent electronic transport 
through an anisotropic magnetic molecule using a simple 
microscopic model that provides spectral properties sim- 
ilar to those of a multi-ion single-molecule magnet but 
also goes beyond the giant spin model. Our results show 
that at very low temperatures certain transitions are sup- 
pressed when an external in-plane magnetic field is tuned 
to a special direction and value. This suppression can be 
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Figure 10: Plot of the differential conductance, dl/dVb, as a 
function of bias voltage at b — 63/2 (dashed line) and for 6 = 
(solid line). The differential conductance is given in arbitrary 
units. 
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Figure 11: Gray plot of the differential conductance, dl/dVb, 
as a function of bias and gate voltage at 6 = 0. The gate 
voltage is shifted by V s = 13.175 V, as explained in the main 
text. 
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Figure 12: Gray plot of the differential conductance, dl/dV 
(in arbitrary units), as a function of bias and gate voltage at 
b = 63/2. The gate voltage is shifted by V s = 13.175 V, as 
explained in the main text. 



interpreted as the result of destructive interference of a 
topological phase - the Berry phase - as the total mag- 
netization of the molecule transits between two opposite 
orientations. 

The fact that we have to rely on very low temper- 
ature (in the submilikelvin range) to visualize features 
attributed to the Berry-phase interference indicates that 
they will be difficult to observe experimentally. However, 
for real single-molecule magnets with a multi-ion core, 
the energy scales controlling these features could easily 
increase by an order of magnitude. Thus, in practice, the 
situation should be more accessible than what our model 
calculations indicate. When searching for Berry-phase ef- 
fects in transport, one should consider studying molecules 



with a large tunnel splitting (which most likely will re- 
quire a relatively strong in-plane anisotropy as compared 
to the uniaxial anisotropy). 

We plan to further explore our microscopic formula- 
tion to consider more complex molecules with multi-ion 
transition-metal cores. We are currently extending our 
approach to investigate Berry-phase effects in the elec- 
tronic transport of mononuclear (single- ion magnets), 
which consist of just one rare-earth ion within a ligand 
cage. These systems have been synthesized with success 



in recent years (see for instance Refs. l5D - fM]) . 
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